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ABSTRACT 


A mesoscale investigation of a series of tornadoes and 
thunderstorms that passed through the NSSL mesonetwork in 
Oklahoma, on 10 June 1967, between 1700 and 2300 CST, was 
somaueted: Utilizing upper air data provided by NSSL, the 
divergence, vertical motion, and energy fields were computed. 

A finite-difference technique for computing and 
smoothing divergence was developed. The vertical motion was 
then computed by means of the kinematic method. Both the 
total derivative and the local rate of change of static 
energy were computed. Using the values of the local rate 
of change, prognostic fields of static energy and an energy 
index were obtained. Finally, forecast energy indexes, 
divergence, and vertical motion fields were compared to the 


observed locations of tornadoes and thunderstorms. 
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I. INTRODUCTION 


Research in meteorology is normally related to one or 


more of three general scales of atmospheric motion. Al- 


though most meteorologists accept the division of the at- 


mosphere into three scales of motion, they do not agree as 


to the exact definitions and dividing lines between these 


scales. 


The Glossary of Meteorology (1959) defined these 


scales as: 


1 Lag 


SyanOpelc Ot Gy Glrolle mMeteOToOLoyy, Which 1chers eGo 
the migratory high and low pressure systems of the 
lower troposphere, with wave lengths of 1000 to 
2500 kilometers; 

MesomecteoroLlogy, which 1S) Concermed wren the deuce 
tion and analysis of the types of major weather 
phenomena that are small enough to remain undetec- 
ted within a normal observation network;! and 


Micrometeorology, which deals with the observations 


and explanations of the smallest-scale physical and 


dynamic occurrences within the atmosphere and is 
normally confined to the surface boundary layer of 


the atmosphere. 


The majority of meteorological research performed has 


been done on the synoptic scale due to the world-wide 


synoptic network of surface and upper air reporting stations. 


1 
Fiedler and Panofsky (1970) defined the wavelengths 
for the mesoscale as lying between 20 and 500 kilometers. 
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In addition, a large amount of research has been undertaken 
on the micrometeorological scale with researchers being able 
to set up dense surface networks to investigate areas even 
smaller than one square kilometer without excessive diffi- 
culty. 

It is in the area of mesometeorology that research has 
lagged. This has been primarily due to the fact that the 
Synoptic data network is not dense enough, even in the 
United States, to allow adequate resolution of such weather 
phenomena as tornadoes and thunderstorms. As a result, 
early mesoscale investigators, such as Fujita (1955) who 
investigated a system of squall lines in the central United 
States, were forced to use synoptic data. Fujita attempted 
to overcome the inadequacy of the synoptic data by converting 
time cross sections to space cross sections. Nevertheless, 
Fujita noted that mesoscale networks were essential if re- 
search in mesometeorology was to advance. 

Due to the efforts of Fujita and others, a surface 
mesonetwork was established in central Oklahoma by the 
NoeLoOnal Severe Stomms Project @(iooP) in the sprimngsot 1900”. 
This network, referred to as one Beta network, was estab- 
lished in order to obtain data on the severe thunderstorms 
and tornadoes which are prevalent in central Oklahoma in the 
spring. Although the network consisted of some 36 reporting 
stations located from 10 to 15 miles apart, it was evident 
that a knowledge of surface conditions was not sufficient. 

Finally, in 1966, the National Severe Storms Laboratory 


(NSSL) of the Environmental Sciences Service Administration 


EO 





(ESSA) augmented the Beta network with a nine station upper- 
air mesonetwork covering south-central Oklahoma. The nine 
stations, located from 40 to 50 miles apart in 1966 and 
1967, took upper-air soundings every 90 ee during 
periods of strong conwmectmve activity in the spring and 
early summer. This resulted in greater in-depth analysis 

of mesoscale phenomena over more extended periods of time 
than previously possible. 

In 1968, NSSL relocated the nine stations such that 
the distance between them was approximately one-half what 
it was previously. This resulted in a more dense upper-air 
reporting network, but at the sacrifice of reducing the 
area available fOr Inves#a cation. In 1970), ESSA» becancea 
department under the National Oceanic and Atmospheric Admin- 
istration (NOAA). 

The inherent errors contained in present day radiosondes 
have led to questions concerning the reliability of the data 
obtained from the more dense network. As a result, some 
researchers have limited themselves to the extensive data 
obtained by NSSL during 1966 and 1967 for detailed studies 
Geesevierc Storms and their associated weather phenomena. 

One such set of data was used for this study. 

On 10 June 1967, convective activity was present through- 
out the NSSL network. As described in Storm Data (1967), 
this convective activity resulted in scattered but frequent 
tornadoes, thunderstorms, hail, high winds, and heavy rains 
during the observation period from 1700 CST through 2300 CST 


on 10 June 1967 (2300 GMT, 10 June 1967 through 0500 GMT, ll 
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June 1967). Indicative of the severity of this weather is 
ewe Ltollowing report contained in Storm Data (1967), for 
the central part of Blaine County, Oklahoma, between 1930 
eoreand 2200 CST, on 10 June 1967: 
"Heavy thunderstorms moving northeastward across the 
County caused heavy wind damage in the Hitchcock- 
Watonga-Greenfield area. One tornado was observed on 
the northwest side of Watonga about 8:30 pm. Three 
Other tummelswaloft were sighted near the city and may 
have possibly passed low over the city during the 
blinding rain. A 1/3 mile-path in the Hillcrest Heights 
addition in northwest Watonga left a demolished hanger 
and five smashed airplanes, 2 homes completely destroyed, 
8 homes with extensive damage, and most others damaged 
to some extent in the city. Roofs of downtown businesses 
were damaged and windows blown out. Three trailer homes 
Wwenbe Gdemolashed. A Second tornado is believed Co havc 
struck the OG&E transformer station 6 miles east of 
Watonga because the 30 tons of steel structure was 
Intted up into the air and dropped 150 feet to the east: 
Nearly every farm homestead in the area received heavy 
outbuilding damage and some lost homes and heavy steel 
buildings." 
Because of the intensity of such storms and the damage caused 
by them, it 1s imperative that research be continued in meso- 
meteorology. 

iInverdéer to mmvesticate arinospher1¢ pienomen ae Get niene 
on the mesoscale, a 12X10 grid was superimposed on the NSSL 
network encompassing an area from 96°45' to 99°30' West and 
mom 55 45" to 36°00 North. in grid sections of 5" fone — 
[iGesand Latitude Tespeetively. —Figure lL shows Ene waverae 
grid with the dashed lines enclosing the reduced grid area 
Wea lized for Most computations. in this study due To theyre 
differencing schemes employed. 

Utilizing the data for the five time frames commencing 


With 2500 GMI, 10) June T9¢7emand spaced 90 minutes apart. 


PISmSpiGvaratt tempted to: 
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Fig. 1. Mesoscale grid network used in study. Dashed 
line encloses area utilized for most computations. Upper 
air reporting stations are indicated by@. 
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Describe the overall convective activity occurring 
throughout the area during the period under inves- 
tigation; 

Develop a method for obtaining realistic vertical 
motion fields; 

Determine what knowledge could be gained from the 
local rate of change of static energy; and 
Correlate the divergence, vertical motion, and 
changes of static energy fields with the areas and 
times of tornadic activity in order to determine 


what mechanisms were present. 
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Beno rNOE Prec oLluemuen Al O000 GMT, 11 June 1967 


To aid in understanding what has occurred on the meso- 
scale, it is helpful to look at the synoptic weaat ion 
during the same time frame. During the period under inves- 
tigation (2300 GMT, 10 June 1967 through 0500 GMT, 11 June 
1967) a synoptic surface and upper air sounding in the NSSL 
network was taken at 0000 GMT, 11 June 1967, at Oklahoma 
City, Oklahoma. This sounding, which took place during the 
first hour of the time period, was taken at a time of marked 
tornado, thunderstorm, and hail occurrence in the northwestern 
Bpection of the mesoscale network. 

Eugunes 2mthrouem 5 represepewtne Synoptic SsitUation 
over the south-central United States at 0000 GMT, 11 June 1967, 
as analyzed by the National Weather Service. Figure Z de- 
tails the surface situation. A stationary front was located 
in northwestern Texas and western Oklahoma, running in a north- 
easterly direction. A low pressure system with a weak wave 
in the front was present in the panhandle area of Texas. The 
wind ahead of the front was from the southeast bringing warm, 
welatively Most alr trom the Gulf vor Wie coe nis camp. 
Sosehved from the temperature and dewpolme. temper acuress plo. 
ted in Figure 2. The dewpoint temperatures behind the front 
were much lower, indicating warm dry air was advected behind 
the front from the deserts of New Mexico and Arizona. Since 
the temperatures on both sides of the front were comparable, 
the primary way to delineate the frontal location was by 


means Gtetwe wind shitt waeressethe front. In this case wemene 
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was a marked veering of the wind from the southeast to west- 
meuthwest, providing a clear indication of the front's loca- 
tion, and also indicating the presence of a low level 
memvergence zone in the vicinity of the font. 

Figure 3 outlines the synoptic situation at 850 mb. The 
front was still evident at this level with southerly flow ad- 
vecting warm, moist air (mixing ratio of 10.0 g/kg) into the 
area ahead of the front with westerly flow advecting warm, 
Gry air (mixing ratio of 5.8 g/kg) behind the front. The 
isotherm pattern showed a definite influx of warm air into 
the mesonetwork area. The isoheight pattern was tightest 
over the center of the mesonetwork and indicated difluence 
in the northern section. Finally, the wind speed was rela- 
tively high over the mesonetwork. | 

Figure 4 shows the height and temperature patterns at 
700 mb for the same area and time. The wind direction was 
from a southwesterly direction throughout the entire area. 
The temperature patterns indicated that warm air with a 
mixing ratio of 6.0 g/kg was advected into the NSSL network. 
This compared with a mixing ratio at the surface of 13.1 
g/kg and at 850 mb of 13.5 Byer although some areas at the 
surface had mixing ratios as high as 15.0 g/kg. There was 
an indication of a convergence zone northeast of the meso- 
network. In addition, the maximum wind speed of 40 knots 
was recorded in the mesonetwork. 

Frgure 5 depicts the aummp hereht and temperature pats 
terns. The wind, as at 700 mb, was predominately from the 


southwest; however, it was generally stronger than at 700 mb. 
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A weak ridge had formed over the central and eastern por- 
tions of Texas and southeastern Oklahoma. The Mixing rales 
throughout the general area was of the order of 0.7 g/kg. 
However, in the vicinity of the mesonetwork the mixing ratio 
was 2.5 g/kg. In addition, a convergence zone was indicated 
in the area of the mesonetwork where the wind speed decreased 
from 40 to 25 knots with a weak divergence zone to the north- 
east of the mesonetwork. Finally, a tongue of cold air was 
moving in toward the mesonetwork from the west. 

The overall conclusion, based on these synoptic weather 
charts, was that in the vicinity of the mesonetwork, cold 
air with relatively high moisture values was located at 500 
mb over warm-air with high moisture values. This, coupled 
tach Glie Fact that convergence was indicated in the vicinity 
of the mesonetwork at all four levels, resulted in a highly 
unstable condition leading to convective overturning. The 
fact that severe weather was present in the mesonetwork at 


this time lends support to these conclusions. 
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Dl. SIN EomeneAl IVE PROCEDURES 


A. DATA ANALYSIS 

Mimonder sto dcecomalishn the objectives of this study, 
it was necessary to analyze three parameters for each of the 
five sounding times (2300, 0030, 0200, 0330, and 0500 GMT) 
and 17 levels (surface to 150 mb at 50 mb intervals). Util- 
izing the computerized data provided by the NSSL mesonetwork 
for 10 June 1967, wind direction, wind speed, and static 
emergy were plotted on the 17 isobaric surfaces for each of 
Bicmmine Stations for ¢ach time period. 

A problem that developed in the plotting and analysis 
of the NSSL data was a combination of either the inability 
of some stations to launch a radiosonde due to severe weather 
conditions or, at some times, the radiosondes failing to 
reach the 150 mb level. This was especially true during the 
last three time periods when Watonga (WAT) failed to launch 
due to tornadoes and thunderstorms in the general vicinity. 
In addition, at 0330 GMT, Tinker (TIK) also failed to launch 
a radiosonde and Cordell's (COR) radiosonde only reached 350 
mb. As a result, the northern section of the mesonetwork 
was almost entirely without reports necessitating subjective 
analysis based on extrapolation from previous times and from 
those reporting stations to the south. Unfortunately, all 
tornado and thunderstorm activity during the period of this 
study occurred north of a line from Cordell to Tinker. 


Therefore, the computed data for the important northern 
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section of the mesonetwork was obtained from extrapolated 
analyses and hence was subject to greater errors than the 
other sections. 

Once all streamlines, isotach, and static energy iat 
terns were analyzed, 120 data points were extracted for each 
level and time frame. In addition, a fourth parameter, sur- 
face pressure, was plotted and analyzed in order to take into 
account the surface effects upon the wind direction and speed. 
Mmiall, 31,200 data points were extracted for the study. 

In the case of all computed parameters, the centered 
@enate difterencing technique resulted in a reduced 10X8 erid; 
Pie parameters listed in this chapter are detailed in the 
computer program CGltained) 1neAppendix A. 

Space and time continuity were maintained throughout as 
consistently as possible. The continuity in space was fa- 
cilitated by analyzing at 50 mb intervals which were chosen 
Because the NSSL computerized data provided the information 
at these levels. Also, this allowed many small-scale motions 
to be picked out and analyzed that otherwise might have been 
missed if a larger vertical increment had been chosen. 

All data was plotted and analyeen as on station data, 
even though in one case the launch times at the nine stations 
varied as much as 25 minutes. On the average the launch times 
varied only by nine minutes. In addition, the balloons ob- 
viously did not ascend in a straight vertical profile. 
Fankhauser (1969), felt that such factors as the lateral 
displacement of the balloon from its launch site, the varia- 


bility of launch times, the differing rates of balloon ascent, 


23 





fad the time lag of the dafferent pressure levels from launch 
time must all be taken into account in performing mesoscale 
analysis. He felt that neglecting these factors would in- 
troduce extreme errors. To compensate for Sheae factors, 
Fankhauser utilized a linear interpolation formula of the 
form 


ceo + (Et | og - x) (1) 


where X is the variable under consideration and the sub- 
menipts | and 2 refer to adjacent soundings in time t.. This 
approach was investigated for use in this study, but it was 
@ecided not to use it for two reasons: 
1. it approximately doubled the amount of data ex- 
traction required; and 
2. It eliminated one time frame from utilization in 


the study. 


Je DIVERGENCE “GOMBUTATIONS 


Velocity divergence, which may be expressed as 


DIVa VY a Oe (2) 
was computed as it is a required parameter for any discussion 
of convective activity on one mesoscale where computation of 
vertical motion by the kinematic method is used. 

A major problem in computing divergence is due to the 
fact that the individual terms are normally of approximately 
equal magnitude but of opposite sign. This results in the 


divergence being one order of magnitude smaller than the in- 


dividialetenns. This Can yestecein large errors 1 either 
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the streamline and/or the isotach patterns are incorrectly 
analyzed, or if there is any error in the original data. 

Although equation (1) was not utilized, it was realized 
that the points raised by Fankhauser were valid, at least to 
some extent. In order to compensate for the deviations pre- 
viously listed, the divergence field was smoothed in both 
eae horizontal and vertical. 

The actual computation of isobaric divergence was done 
utilizing a centered, finite-difference, trapezoidal area- 
averaging technique based on Gauss! Theorem. It may be ea 
pressed as follows: 

oo 1 

VeVo = - 1 Oc 61 (3) 

a ) 

where V°V>2 is the area-averaged isobaric divergence, A is the 
area computed using the trapezoidal rule, Ch 1s the normal 
component of the wind at the boundary, and 61 is the incre- 
mental length along each side of the boundary. Figure 6 de- 
picts a representation of the grid used. Area averaging in 
the horizontal succeeded in smoothing out many of the irreg- 
ularities normally present in elementary divergence computa- 
EONS . 

Once the divergence field had been smoothed in the hor- 
izontal, it was necessary to determine whether or not to 
smooth the divergence in the vertical. Computing divergence 
through all 17 levels with no vertical smoothing resulted in 
patterns as shown in Figure 7. In the case presented, the 
divergence pattern was not vertically consistent in that there 


were adjacent extremes of both convergence and divergence. 
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aguithary anterior point (1,7). 
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Pesed on this, plus the £act that the divergence at any 
point must also be affected by the divergence at surrounding 
points, it was felt that a vertical smoothing should be per- 
formed and compared to Figure 7 in order to determine if a 
more realistic divergence pattern could be obtained. 

Numerous vertical smoothing techniques were tried and 
discarded before a satisfactory one was finally obtained. 
The resultant vertical smoothing technique employed a three 
dimensional, space-centered, finite-difference smoothing 
scheme as depicted in Figure 8. 

In the actual computation, a weighting factor was as- 
Signed to each point with the heaviest weight being given 
Moewene point being smoothed. Since the divergence field had 
already been smoothed in the horizontal, the next heaviest 
weight was equally given to the points immediately above and 
below the poamt to be smoothed. Finally, the horizontal 
points on each side were assigned the smallest weight. Var- 
lous combinations of weights were tried to determine which 
combination would result in the most representative and 
Realistic divergence patterns. fhe weighting factors fine, 
chosen were those given in ieee Si: 

The smoothing technique was modified at K=1 (950 mb) 
and K=17 (150 mb) as well as at the lateral borders since a 
Space-centered smoothing scheme was not possible. This re- 
sulted in the elimination of one weighting factor of 3 at 
K=1 and K=17, and one weighting factor of 1 on each of the 
lateral boundaries. The results of this smoothing of the 


divergence field are shown in Figure 9 which depicts the same 
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Fig. 8. Example of three-dimensional 
smoothing technique utilized in smoothing 
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grid point as shown in Figure 7. Comparing these two fig- 
ures, it was obvious that the extremes in the divergence 
field of Figure 7 had been damped in Figure 9, while the 
vertical motion patterns in both figures remained basically 
the same. It was therefore decided that the vertical smooth- 
ing of the divergence field was useful. 

A final computation involving the divergence was to sum 
the divergence in each vertical column to see if the sum of 
the divergence was small. No values approached zero, but 


4 


were of the order of 10 *’ to 10 ° per second. If one used 


these values to compute the vertically averaged divergence, 


> see) ali) © per “Sec= 


the results would be of the order of 10 
ema. These values were deemed acceptable since an artificial 
med, ague to lack of data, had been placed on the atmosphere 
at 150 mb while convection may well have extended above this 


level. 


e. VERTICAL MOTION COMPUTATIONS 

Several vertical motion schemes have been developed for 
computing values of vertical motion (omega) on the synoptic 
scale. Whether or not these schemes would work on the meso- 
scale was the question. Sines constant pressure surfaces 
were utilized in the analysis, this automatically eliminated 
computation of vertical motion by any of the isentropic 
methods. The thermodynamic methods were eliminated because 
of the difficulty in obtaining values for diabatic heating. 
The quasi-geostrophic form of the omega equation was not used 


because it required a Rossby number of less than 0.1, and on 


Sa 





the mesoscale it is not uncommon to obtain Rossby numbers 
emeater than 2.0 in areas of strong convective activity. 
hrtcreclininaeioneor the abeve methods, only the kine- 
matic and vorticity methods for computing vere motions 
were left. The vorticity method, as discussed by Coleman 
moo9), was eliminated due to the result that computed 
values for vertical motion were not within the realm of 
acceptability on the mesoscale; that is, the values obtained 
became so large that they could no longer be believed. 
Mitts Dy ene process Of Clifination, the kinematic 
method of computing vertical motion was chosen. However, 


eae Kinematic method involves integration of the continuity 


equation: 
Ue oo (0 (4) 
ap 

which contains many disadvantages. These disadvantages may 


be categorized as follows: 


fey Etetenr sane upper On lower se Mld anya come leGionle lls aaecm 
Giiring ed. 
7 1G as computed tromecdavencences Wiel. Ele Gis acme 


tages discussed in the’previous section; 


3. It is very poor in regions of ill-defined winds; 
and 
4. It tends to yield excessively large values, even 


in areas of good wind data. 
Before the kinematic method could be used to compute the 
vertical motion it was necessary to provide solutions to the 


disadvantages listed above. An assumption of w=0 at 150 mb 
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was chosen in order to allow computation of a linear error 
term. This appeared to be an adequate assumption based on 
an inspection of the analyses which showed that the tropo- 
pause was located in the vicinity of 150 mb and at wivemitadc — 
popause vertical motion is generally assumed to approach 
zero. It was realized, however, that some vertical motion 
would probably penetrate the tropopause due to the strong 
convective activity that is known to be associated with 
Severe thunderstorms. But since insufficient data (only 
mye or three out of nine stations) was available above 150 
mb to allow for a realistic analysis, this assumption was 
considered necessary. 

The lower boundary condition of w=0 at the surface was 
not utilized because of the large amounts of convective ac- 
tivity and vertical motion that would normally be anticipated 
in an area of tornadoes and thunderstorms. Instead, the 
vertical motion at the surface was computed utilizing the 
terrain effects of the u and v components of the wind and the 


surface pressure, that is: 


— op op = ie 
"sfc Per eae oy y “Ps fe (5) 


where VeVi Poe. is expressed in cartesian coordinates. 

The second and third disadvantages were minimized in 
this case due to the averaging techniques employed for the 
divergence in both the horizontal and vertical plus the fact 
that there were nine upper air reporting stations in the 
mesonetwork. 

The last disadvantage presented the greatest degree of 
difficulty and the chief reason why most meteorologists have 
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not used the kinematic method in the past. Gruber and 
O'Brien (1968), and O'Brien (1970), have investigated this 
peoblem in great detail. O'Brien (1970) found that case 
studies for intense mesoscale convection dihengeraced the 
reliability and usefulness of an adjustment scheme for the 
kinematic method which was based on the premise that the 
errors in the divergence field were a linear function of 
pressure. 

Sete peeawithiethesmimtesgrated form of the continu ty 
equation: . 
oe (6) 


WWWV4 
where the integral average is approximated by veV, (the 


Mertical average of the divergence field between two levels) 


meme Values of omega were successively computed, that is: 


9 + DIV. 


f£ 
20 . C | (Pip. - 900). (7a) 


4 


fea 


Utilizing the value for omega at 900 mb to obtain the value 


Of omega at 850 mb, yields: 


DIVore, + DIV 
: 850 900 
“850 ~ 900 * | se | (50) 
SEE so0 8 2 ote | @ ae 3 00) 
ae 2 
oe PIV 50 + DIVo 99 (7b) 


Continuing this cumulative procedure through to 150 mb gives: 


+ DIV. 


900 1 es 
5 | (ees - 900) 


+ 


DIV 
Ww = Wy 
150 SEC 


= 25 (DIV) ue DULY! He eV. 


00 £50 
34 


(850+200)? (7) 





where Ge 200) means the sum of all of the divergence 
terms between 850 mb and 200 mb. Now, it was assumed earlier 


that Wicg 0, thus, equation (7c) may be rewritten as: 


DIV ae WL 
: 900 sfc : 
ee oaerne i z (Eg 900) 
+ 25 (DiVo99 * DIVi sg + 2EDIV g59.209))- (7d) 


Dividing equation (7d) by 25 mb in order to place the equa- 


f?on in units per second yields: 





0 = “sfc A rego @!’ ste (P - 900) 
Z5 50 Se 
og pe iso) 2 DEY casos zn) . (7e) 


tse Wisg TCC alavenwias 2 C110) ee © all divergence CennSe Wetec 
exact, and if there was no computer round-off error, then 
eaimation (€7e) would be identically satisfied. _It 1s obvious 
that this is not the case and thus (7e) is only approximately 
satisfied. The difference between zero and the value ob- 
Gimed by the above relatvonship 15 the error term. Calling 
this divergence error term DVERT, and setting it equal to 


the previous relationship yields: 





m Dt ayay eee Ln 
Wy ke 900 Sie : 
pvert = -§£2 4 |— 900 sfc | cp .. - 900) 
eC) Be ee aso AU) . (7£) 


Now dividing DVERT by the sum of the multipliers of the 
divergence terms and the vertical motion at the surface, or 


Salada deme by, (P.¢.7150)x0.04, resulted in an average error 
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for the divergence throughout the vertical. This corresponded 
to O'Brien's (1970) statement that the errors in the diver- 
Menace t1ela were a linear function of pressure. Subtracting 
this error from the computed averaged divergence at each 

level resulted in a corrected divergence which was then used 
in the computation of the vertical motion. The results of 
M@tilizing this method can be seen from Figure 7 and Figure 9. 
In addition, the average values of omega at 150 mb in this 
study were 10.’ mb sec.” or 8.6x10 ° mb day *, well within 

mie acceptable tolerances established prior to beginning 


this study. 


D. Slat TC ENERGY COMPUTATIONS 
Total energy may be defined as 


V2 
Ee ee (8) 


where CS 1 the specitic heat of alr ae ConStanepnessume. 
T the temperature, gZ the geopotential, L the latent heat, 
q the specific humidity, and V the scalar velocity. The 

substantial derivative of total energy is given by file. shee 


lationship: 


2 
we (CT + e2+lg+ 5) = Gira he E (9) 


where dh/dt represents the sum of sensible heat, evaporation, 
amd radiation, @ 1S the specitic yolume, dp/dt 15 the Vocal 
pressure change, and F represents frictional heating. Now, 
assuming that the sum of the terms on the right of equation 


(9) are sufficiently small, one may write 
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d C T+ +) Lacs Me = 0 
at {C, gz aa) = 2 - (10) 


Since static energy is defined as 


E. = c,! cee LC (11) 


equation (10) may be written 


d d Me 
Since d/dt(V*/2) is a relatively small term, as verified 
by utilizing observed winds, the static energy was assumed 
to be conservative. | 


The conservative nature of E. may be further understood 


by writing 


tT 


0 
4, pores + 7. + lL 13 
p ba & q) (es) 


C 
p 


=| — 
| 
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and recalling that 6, is conserved for both dry- and saturated- 


E 
adiabatic displacement. It follows that the difference in 
static energy between two levels is directly related to the 
BOrential instability of the layer. 

In addition, the normal instrument tolerances in a radi- 
osonde allow the additional assumption of q = w, where w is 


the mixing ratio. Thus, the static energy may be expressed 


by the following approximate relation: 


E. x oH aero edly (14) 


The first item investigated was the total derivative of 
Seavicecnergy. Utilzingeamspace- “and time-centered finite 
differencing scheme, the total derivative of static energy 


was obtained using the expansion: 


oi 
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ices sox 8 cy - op 


JE 
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The results obtained from computing the total derivative 


were of the order of 10. 


‘ m2(s?-1.Shr) ~~ 


The total deriva- 


tive of static energy was computed to check how well the 


finite differencing scheme employed conserved static energy. 


Assuming dE _/dt=0, thesleeal rate 
energy was computed from the advective 


mare of change of static energy can be 


The values obtained for the local 


JE JE 


ae fo) ES) 
Sot wo ox dy iia 


JE 
aac 
dp 


Static energy using equation (16) were 
magnitude as those values obtained for 


@2) Static Energy. 


OL -Ghamge= OF —S tate 
Cems nem loea. 


expressed by: 


JE 


+) —— 


(16) 
dp | 


V9 
S 


fate. Of eGhiqmecmon 
of the same order of 


the total derivative 


These values were then used as a prognos- 


tic tool by taking a single forward time step to obtain the 


matic energy at time t + AL, where At 


may be expressed in equation form by: 


Pteat ~ Pe * iat 


No boundary conditions were required in computing E 


t| At 


TS ke So nouns. nes 


(17) 


CPE 


as the forecast was confined to the interior grid points. 


It was realized that the prognostic values of E 


Parsee would 


not exactly equal the actual values due to the assumptions 


previously made. 


nostic values could be used to predict 


However, it was desired to see if the prog- 


an accurate energy 


index based on the premise that the generated errors of the 
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prognostic static energy would cancel each other out in the 
computation of the energy index. 


The Energy Index was defined by Darkow (1968) as 


ya ee (Ercog - Eres) (18) 


with values in joules om indicating the following: 

Pee tom 4 62 thunderstorms possible; 

Pema --  1SOlated Severe thunderstommede— 

tivity possible; and 

3. <-8.4 -- severe thunderstorms and associated tor- 

Mado activity are highly probable. 
The above criteria are conditional upon an adequate trig- 
Bering mechanism being present or forecast to release the 
mocential instability. 

The prognostic Energy Indexes obtained by subtracting 
mae LOLeCast Static energy at 850 mb from that forecast for 
muOo mb, were compared to the actual Energy Indexes to deter- 
mime it they were in general agreement. Finally, the static 
energy computations were compared with the divergence and 
meTeical MOtion patterns to determine 1f any Cornelatiron 
existed. Results of these comparisons are included in the 


following section. 
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IV. RESULTS 


Ij iiciimibatcmtmcmamselussion of the results ofuthis 
Study., each time frame is looked at separately and in se- 
amence beginning with 2300 GMT, 10 June 1967. All figures 
in this chapter are on a 12X10 grid, with parameters calcu- 
lated on the 10X8 grid extrapolated one grid distance in 
each direction on the boundaries. This expansion of the 
reduced grid was necessitated by the fact that many of the 
tornadoes and thunderstorms occurred on the northern bound- 
fey OL ethe 12X10 grid, ,and utilization of the 10X8 grid would 
Besuit in these features being excluded. 

The space-centered finite differencing schemes used 
throughout this study along with the original data and anal- 
ysis problems resulted in the border areas being the areas 
where the greatest errors occurred. To prevent any bias 
from being introduced into the analysis, it was decided that 
the locations and times of the tornadoes and thunderstorms 
Should not be known until all analysis and data extraction 
had been performed. As a result, the areas of greatest 
tornado and thunderstorm activity were also the areas where 
the greatest probability of errors existed. 

The 750 mb level was used throughout this section to 
represent the divergence and vertical motion patterns for 
all time frames. This level was chosen because it was il- 
lustrative of the low level air involved in the convection. 

All tornado and thunderstorm locations indicated on the 
figures in this chapter were obtained from either Storm Data 
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(1967) or zero degree elevation radar pictures taken by the 


WSR-57 radar located at Norman, Oklahoma. 


A. THE MESOSCALE SITUATION AT 2300 GMT, 10 June 1967 

‘Figure 10 shows the actual Energy Index for the first 
time frame. From the values shown in this figure, the en- 
tire mesonetwork was an area of probable tornado activity 
with all values more negative than -12, and with a rela- 
tive maximum probability located in the northwest quadrant 
where a tornado did occur at 2309 GMT. A thunderstorm. was 
miso contained in the area of maximum probability. 

Figure ll shows the divergence and vertical motion 
patterns at 750 mb for the same time. The tornado and thun- 
derstorm were seen to occur in an area of ascending vertical 
motion and strong horizontal convergence. In addition, 
there were relative maximums of both convergence and ascend- 
ing vertical motion ahead of the tornado and thunderstorm. 
Behind these features, there were areas of divergence and 
ilescending vertical motions. Both oi these results éon-— 
formed to models previously developed for tornadoes and thun- 
derstorms (cf. Browning and Ludlam, 1962). It may also be 
noted from Figure 11 that it was not possible to exactly 
locate the tornado based on the divergence and vertical mo- 
tion patterns as a tornado is normally considered to be a 


sub-mesoscale phenomena. 


B. THE MESOSCALE St iVATTeGNeEOR 0050>6GMI, 11 June 1967 
Figure 12 shows both the actual and the forecast Energy 


Indexes for 0030 GMT, 11 June 1967. The actual Energy Index 


4] 








ee 


fe, LO. Energy Index for 2300 GMI, L0°-June 1967. Alt 
Values are in joules gm~*. w™indicates tornado activi cy, 
q indicates maximum radar return, and ©) indicates 
secondary radar return. 
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foo, Jil. Divergence (ene aioe sec and Vemeical motrvon 
(——)x10 ° mb sec * at 750 mb for 2300 GMT, 10 June 1967. 
V indicates tornado activity, g@ indicates maximum radar 


return, and Q@)) indicates secondary radar return. 
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Fig. 12. Actual (——) and forecast (---) Energy Indexes 
for 0030 GMT, 11 June 1967. 
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indicated once again that the entire mesonetwork was a prob- 
able tornado and thunderstorm activity area. There was a 
relative maximum area of probability in the northwest and 
a relative minimum in the northeast. 

iMiemprognostic Enérgy Index also indicated that the 
entire mesonetwork, with the exception of the western sec- 
tion, was a probable tornado and thunderstorm area. An 
area of relative maximum probability was located in the 
northwest with the area of minimum probability situated in 
mre west. A seconde minimum was located in the northeast. 

The two areas of relative maximum probability compared 
favorably as did the actual Energy Index minimum with the 
forecast Energy Index secondary minimum. The tornadoes and 
mrondaerstorms both occurred in the relative maximum areas 
Of Probability, thet 15, in areas smaller than -loyjfowlles 
gm 

There was generally good correlation between the sewo 
indexes, although the western area of the forecast Energy 
Index was an exception. Figure 13, the divergence and ver- 
tical motion patterns for the same time period, showed that 
the area of minimum probability in the forecast Energy Index 
corresponded to an area of Sstrome descending vertical netion 
and divergence. The secondary area of minimum probability 
in the forecast Energy Index also corresponded to an area of 
divergence and weak descending vertical motion. 

Both Energy Indexes correlated well with the vertical 
motion patterns in the vicinity of the tornadoes and thunder- 


storms. However, there appeared to be a stronger relationship 
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between the prognostic Energy Index and the divergence and 
Mertical motiom patterns than between the actual Energy 
Index and the divergence and vertical motion patterns. This 
was not considered unusual when the method of obtaining the 
prognostic Energy Index was reviewed. The prognostic Ener- 
gy Index for t+At was obtained by adding the static energy 
memtame t to the local rate of change of static energy at 
time t; and the local rate of change was heavily dependent 
mr the u, v, and w components of the wind. 

From Figure 13 it can be seen that the tornadoes ane 
epUunaderstorms once again were located in areas of ascending 
vertical motion and convergence. Ahead of the tornadoes 
and thunderstorms there was ascending vertical motion and 
convergence, while behind them there was descending vertical 
motion and weak convergence changing to divergence. Once 
again, it was not possible to determine the exact location 
Gf the tornadoes based on the divergence and vertical mo- 


Peon patterns. 


er: THE MESOSCALE SITUATION FOR 0200 GMT, 11 June 1967 

The observed Energy Index and the forecast Energy Index 
for 0200 GMT, 11 June 1967, are contained in Figure 14. The 
observed Energy Index indicated that the whole mesonetwork 
area was a probable tornado and thunderstorm area. A rela- 
tive maximum area of probability was located in the northwest 
sector with a relative minimum area of probability located in 
the northeast section. The tornadoes and the thunderstorms 


actually occurred in the Morenmepetween the two extremes. 
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Fig. 14. Actual (——) and forecast (---) Een Indexes 
mor OCOUNGME, Li June 1907; 
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The prognostic Energy Index also indicated that the 
entire mesonetwork was a probable area of tornado and 
thunderstorm activity, with the exception of the southern 
mmea, ihe forecast Energy Index indicated a relative maxi- 
mum area of probability in the center of the network with 
two areas of minimum probability, one in the southern sec- 
tion of the grid and the other located in the northeastern 
meerion. As in the Case with the actual Energy Index, the 
tornado and thunderstorms occurred in a highly probable 
tornado and DRUMS Lan area. 

Figure 15 shows the divergence and vertical motion pat- 
terns for the same time frame. The tornado and thunderstorms 
occurred in an area of ascending vertical motion and conver- 
memece. strong ascending vertical motion and convergence were 
ahead Be en although a pocket of divergence and descending 
mentical motion appeared to be in the area. The word ‘“ap- 
peared" was chosen because it was not known what actually 
existed in the area as Watonga did not launch a radiosonde. 
Thus, all data north of Cordell and Tinker had to be extrap- 
olated both in time and space as previously discussed. This 
is obviously an unsatisfactory procedure in dealing with a 


mesoscale investigation of tornadoes and thunderstorms. 


Dd. THESMESO SCALE ST RUALION FORSUSS0 tii 1 June nso? 
Figure 16 outlines the actual and forecast Energy In- 

dexes for the fourth time frame. The actual Energy Index 

indicated that the entire mesonetwork was a probable tornado 


and thunderstorm activity area, with a relative minimum of 
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Fig. 16. Actual (——) and forecast (---)}) Energy Indexes 
hOre0S50° GMI, Ll June 1967. | 
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REopabtiaty located an the cénter of the network and in- 
creasing outward. 

The prognostic Energy Index also indicated that the 
entire mesonetwork was a probable tornado and thunderstorm 
merivieyearea. Inere was a relative minimum located in the 
Same general area as for the actual Energy Index. Two areas 
Piemaximum probability also existed; one in the northwest 
meeinosother an the east-central section of the network. 
The tornadoes and thunderstorms occurred in areas of maximum 
mrobabrlity for both forecast and actual Energy mdeceer 

MicmcOrrelatron Detween tne actual Energy Index andetne 
merecast Energy Index was poor. However, there was a stronger 
Zamacronship between the rorecast Energy Index and the ver- 
imeed! MOtiON pattern, Figure 17, in the vicinity of the €or- 
macgees and thunderstorms than between the actual Energy Index 
a@e ene Vereical motion patterns. 

The tornadoes all occurred in areas of ascending verti- 
Ga motzon and low-level horizontal convergence, with tne 
same features located ahead of them in their projected path 
of movement. Some portions of the thunderstorms, however, 
Heclr rca Gnedredas OF descending Vertical MOuLom ana stnome 
horizontal convergence; an apparent contradiction. 

Re-analysis of the streamline and isotach patterns for 
this time period failed to improve the vertical motion pat- 
terns at low levels. The probable cause for this failure 
was the fact that both Watonga (WAT) and Tinker (TIK) failed 


to launch radiosondes, Watonga for the second successive time. 
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Fig. 17. Divergence (---)x10 ~ sec and vertical motion 
(——)x10-° mb sec™’* at 750 mb for 0330 GMT, 11 June 1967. 
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As a result, it was necessary to extrapolate the entire 
northern section of the network for two time frames as well 
as extrapolate in space from the southern and central sec- 
mrons Of the network. A single time frame extrapolation at 
0200 GMT resulted in an area of probable error as previously 
discussed. Thus, with the large amount of convective ac- 
tivity present during this time frame as evidenced by the 
many thunderstorms and tornadoes, it was not possible to 
adequately portray the actual divergence and vertical mo- 
mone patterns present. All of this points out the fae that 
any numerical investigation of mesoscale phenomena is severely 


hampered when space and time extrapolations are required. 


ES THE MESOSCALE SITUATION FOR 0500 GMT, 11 June 1967 

The actual and forecast Energy Indexes for the last time 
frame investigated are shown in Figure 18. Once again, the 
actual Energy Index showed the entire mesonetwork as a prob- 
able tornado and thunderstorm activity area with the north- 
western section as a relative maximum area of probability 
and the RORenees cern area a relative minimum. 

The prognostic Energy Index contained a maximum area of 
probability in the western section of the grid and an area 
Orerelative Minimum probabilaty In the morth-centrall Ssectrom: 
Both the actual and forecast Energy Indexes had the tornadoes 
and thunderstorms occurring in areas of probable activity. 

A close look at the divergence and vertical motion 
fields, Figure 19, showed many apparent inconsistencies, such 


as strong ascending vertical motion occurring along with 
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Fig. 18. Actual (——) and forecast (---) Energy Indexes 
hore OSGeo bi June 1967. | 
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Mio, 19. Divergence (ee alos sec. and vertical motion 


(—~)x10 * mb sec-* at 750 mb for 0500 GMT, 11 June 1967. 
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Parone miorizamial divergence in the eastern section of the 
@rid,. Although this may appear to be a contradiction, it 
must be remembered that only the 750 mb level is illustrated. 
For this case, there was very strong horizontal convergence 
mom the surirace up to 800 mb, wnere it changed to horizontal 
divergence. The values of the ascending vertical motion had 
already started to decrease by 750 mb and finally became 
descending vertical motion by 650 mb. Thus, no inconsis- 
tencies actually existed. 

Poulin eunhesnrevious two) time periods, Watonga (WAT) 
miebed to launch a radiosonde. In addition, Tinker’s (TIK) 
radiosonde data ended at 500 mb and the data was so dif- 
menent from the dawanaesOenenr —sitatlons sds) CO, De Nigminesius. 
pect. Once again, the entire northern s0Geon analysis 
depended upon a time extrapolation (this time covering three 
time periods) along with a space extrapolation from stations 
Goethe South. 

This was obviously an unsatisfactory method of deter- 
mining the divergence and vertical motion patterns, not only 
for this time frame, but for the previous two time frames as 
well. However, with no data available in the northern sec- 
Pom Or the grid durimeg the Past threestimes, the enoree mwas 
either perform extrapolations or not use the last three time 
Bemaods at all. “Thus, although the results tor these last 
inieee eine periods may be poor compared to what actually oc- 


Cinred wthey are still tar superior to no results sat all. 
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V. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY 

On 10 June 1967, a series of tornadoes and thunder- 
storms passed through the NSSL mesoscale network in south- 
central Oklahoma, causing wide-spread damage and destruction. 
Utilizing the data provided by NSSL, this study investigated 
Meee divergence, vertical motion, and static energy relation- 
Smaps that existed in the NSSL network on that date during 
meee hour period from 23500 to 0500 GMT. 

Divergence was computed using a trapezoidal, area- 
averaged, smoothing technique based on Gauss' Theorem. In- 
westifpation revealed that vertical smoothing of the diver- 
gence field was also desirable in keeping with Fankhauser's 
(1969) hypothesis that for mesoscale data, some form of 
smoothing is necessary. A five-point weighted smoother was 
Utilized to compute the divergence field. 

From the assumption that vertical motion at 150 mb was 
zero, a correction term was obtained that was linearly sub- 
PereEcad wi LOl tie d1iVeroeneemiicldsdt cach point in accordance 
with O'Brien's (1970) hypothesis that the errors in the diver- 
gence field are a linear function of pressure. Utilizing the 
corrected divergence field, vertical motion was computed by 
means of the kinematic method. The resulting computed ver- 
tical motion field was well constrained. 

Investigation of the energy field centered about static 


enerey. Operating on the pramepple of ConServation of energy, 
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the total derivative of static energy was computed using a 
Space and time centered fimite difference scheme. This was 
done as a check on the data and the computational scheme. 
Assuming that the total derivative of static energy was 
zero, the local rate of change of static energy was computed 
Dy setting it equal to the negative of the advective terms. 
Static energy fields were then forecast using the advective 
mianeges and a forward time steep. The forecast static energy 
was then used to compute a prognostic Energy Index by sub- 
@raeting the values at 850 mb from the values at 500 mb . 
iicapmociestle Enenoy Index was compared to Lhe actual 
Energy Index to determine if a correlation existed between 
ie two indexes and between the indexes and the location of 
the tornadoes and thunderstorms. Finally, the divergence 
mca vertical motion fields were compared to the location of 
the tornadoes and thunderstorms and to the Energy Indexes for 


possible correlation. 


JBee CONCLUSIONS 

Based “upon the results contained in’ Section iV, 7the 

following principle conclusions were made: 

i They computationalemethod® develioned tTomsmooth the 
divergence field reswited in a more realistic 
davemecnecemod etc tiie 

2. The corrected divergence constrained the kinematic 
method used in computing vertical motion and re- 


Sulted wen realistie vertireal motion frelds. 
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Lack of data in the northern section of the meso- 
network during@the last three time frames resulted 
in rather poor divergence, vertical motion, and 
Energy Indes patterns fer those times. 

ihe urand vy conponents of the wind and the static 
energy fields should be smoothed prior to comput- 
ing the total derivative and the local rate of 


change of static energy. 


of RECOMMENDATIONS FOR FUTURE RESEARCH 


fine bines with the results and conelusions obtained ain 


this study, the following recommendations were made concern- 


ing further research in the NSSL mesoscale network: 


i 


That a numerical program be developed to smooth the 
u and v components of the wind in such a manner 
that the divergence obtained by using smoothed 
values of u and v is the same as the divergence 
computed using Gauss' Theorem. 


{That a comparison of divergence and vertrcaly moeren 


fields be made utilizing the method contained in 


this study to the method suggested by Fankhauser. 

That burther resecaren be eCOnmducted Into Ute liza 

static energy as a prognostic tool to localize the 
dbeasOn LULULe Probable tornado  acrivi ty. 

That a numerical program be developed specifically 
for the NSSL mesoscale network that would analyze 

data for all desired parameters in order to remove 


Ehe subjective Crromeerom the analysis. In addition, 
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ME CouUla@mciUppLement the objective analysis tech- 
niques developed by NSSL which deal mainly with 
the 500 mb region as described by Inman (1970). 
These are but four of many possible reconmendathone 
for future research in the field of mesometeorology. There 


are still many questions to be answered. 
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APPENDIX A: COMPUTER PROGRAM 
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